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Brazil 

We consider S'-matrix correlation fnnctions for a chaotic cavity having M open chan¬ 
nels, in the absence of time-reversal invariance. Relying on a semiclassical approx¬ 
imation, we compnte the average over E of the qnantities Tr[S'f(h^ - e)S(E + e)]"-, 
for general positive integer n. Our result is an inhnite series in e, whose coefficients 
are rational functions of M. From this we extract moments of the time delay matrix 
Q = -ihS^dSIdE, and check that the hrst 8 of them agree with the random matrix 
theory prediction from our previous paper [M. Novaes, submitted]. 
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I. INTRODUCTION 


Quantum scattering processes at energy E can be described by the scattering matrix 
S{E), which transforms incoming wavefunctions into outgoing wavefunctions. This matrix 
is necessarily unitary, in order to enforce conservation of probability and, consequently, 
conservation of charge. We consider a scattering region (‘cavity’) inside of which the classical 
dynamics is strongly chaotic, connected to the outside world by small, perfectly transparent, 
openings. We assume M open channels, so that S is M-dimensional. We also assume there 
is a well dehned classical decay rate F, such that the total probability of a particle to be 
found inside the cavity decays exponentially in time as The quantity td = l/F is called 
the classical ‘dwell time’. 

The study of time delay requires the energy dependence of the S matrix, as the Wigner- 
Smith time delay matrixii^ is dehned as Q = -ihS^dSjdE. Properties of Q may also be 
derived starting from correlation functions 


Cn{e,M) 



( 1 ) 


where (■) denotes an average over E. For example, the average value of time delay moments 
Mm = ]gTr((5’") can be obtained as^ 
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In the semiclassical regime (when h 0 and the electron wavelength is much smaller than 
the cavity size), we may use the semiclassical approximation, in which elements of the S ma¬ 
trix are written as sums over classical scattering trajectories^. Calculation of energy-averaged 
transport statistics (like condutance, shot-noise, etc) then require so-called action correla¬ 
tions, sets of trajectories having the same total action, leading to constructive interference. 
Using only identical trajectories and ergodicity arguments^'- one can recover some semiclas¬ 
sical large-M asymptotics. Quantum corrections, important at hnite M, can be related to 
non-identical trajectories having close encounters^, and may be obtained systematically-^—. 

The semiclassical approach can also be used to study time delay. Interestingly, in this case 
one can use the periodic orbitsi^ that live in the fractal chaotic saddle of the system-^ (some¬ 
times called ‘the repeller’). This approach was followed It is actually equivalenti^. to 

the one based on scattering trajectories^^, which Berkolaiko and Kuipers used to treat o 
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semiclassically, initially in the large-M limits and later up to the hrst hnite-M corrections^^ 
Another semiclassical approach to time delay, that avoids correlation functions, has recently 
been introducedJ^ 

We hereby advance the semiclassical approach by deriving a formula for correlation func¬ 
tions Cn{€,M) which is a Taylor series in e, with coefficients that are rational functions of M. 
These coefficients are expressed as hnite sums involving characters of the symmetric group 
and Stirling numbers. Our method is an extension of a recently introduced semiclassical 
matrix model for transport statistics^2i 

Statistical properties of Q can also be calculated using random matrix theory (RMT) ^2.^^ 
In our previous paper^ we obtained the average value of general polynomial functions of Q. 
Equivalence between the semiclassical and RMT approach has long been conjectured, and 
was previously known to be true to leading orders in 1/M. The structure of our formula for 
(7^(6, M) suggests that this equivalence holds exactly in M for all polynomial functions of 
Q. We are able to verify this in many cases, but come short of showing it in full generality. 

We remark that the semiclassical approximation provides the energy-dependent correla¬ 
tion functions, which have more information than the energy-independent RMT statistics 
obtained inM. For instance, correlation functions are required in order to develop a semiclas¬ 
sical treatment of Andreev systems ^222 Also, the semiclassical approximation is in principle 
able to go beyond RMT by including Ehrenfest time effects^22“22 These developments are 
outside the scope of the present work, but we hope they will attract attention in the future. 

A last remark about our semiclassical calculation. It is based on an integral over N- 
dimensional complex matrices, and requires that we take the limit 0. This limit is 

needed to enforce that our semiclassical expansions do not contain periodic orbits. It is 
easily taken in the perturbative framework (see Section III.B), i.e. order by order in 1/M. 
However, we cannot rigorously justify it for the exact calculation. The same issue exists for 
transport statistics^ We believe the nature of this limit is an interesting open problem that 
deserves further study. 

This paper is organized as follows. In the next Section we present and discuss our 
results. In Section 4 we present our calculations. They rely on some well known facts about 
symmetric functions and the permutation group, which we have reviewed in our previous 

paper 
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II. RESULTS 


We develop a new formulation for the semiclassical approach to time delay, based on our 
previous work on transport statistics2i. This requires a matrix integral which is designed to 
have the correct diagrammatic expansion, so that it mimics the semiclassical approximation 
to the correlation functions Cn{€,M). 

Solving exactly that matrix integral, we arrive at a formula for Cn{^,M) in the form of 
a Taylor series in e. Let Xa(/^) be the characters of the irreducible representations of the 
permutation group and dx = xa( 1”) be the dimension of such a representation (we have 
reviewed these concepts in our previous paper^). Our formula is 

1 


( T\j\ " ^ V' V' ^ ^ r \ r 


(3) 


where is some complicated function for which we have an explicit form (see Section 
4.4.1), and 

dO dA) 

(4) 


dA) £{X) 

[my =n [M]x =n [^+^ - i]a. 

2=1 2=1 

are generalizations of the rising and falling factorials. 


For the simplest correlation function, explicit calculations suggest that the following 
expression holds: 

°° { Mo An n-1 1 

(5) 




We can establish that the leading order in e is given, for any n, by Cn{e, M) = l + me + 0(e^). 
To leading orders in 1/M we hnd, for example, that 

1 e2 e2(l + 12ze-8e2) 


01 = 


1-ie M2(l-ie)5 
which is indeed in agreement with the hrst 3 orders as computed from ([5]). 
correlation function we do not have a simple formula, but we can show that 

(l-2ze-2e2) e2(4 + 8ze - - 2ie3) 


+ 0(1/M6), (6) 

. For the second 


02 = 


+ 0(1/M^). 


(7) 


(l-ie)4 M2(l-ie)8 

This generalizes some results that appear in the Appendix oS. 

The average value of moments Aim have been computed semiclassically up to the hrst 
few orders in IjMM Using our new exact semiclassical expression for we could 

compute them in closed form as rational functions of M up to m = 8, and check that they 
agree with RMT predictions from!^^. Unfortunately, we could not establish this agreement 
in general, because of the complicated nature of the function 
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III. SEMICLASSICAL APPROACH TO CORRELATION FUNCTIONS 
A. Semiclassical Approximation 


In the semiclassical limit /i ^ 0, M oo, the element Soi of the S matrix may be 
approximated by a sum over trajectories 7 starting at channel i and ending at channel 


^ y (8) 

V ^ H 

The phase is the action of 7 , while is related to its stability. The prefactor contains 
the so-called Heisenberg time, = Mto- 

Consider the correlation function Cn{e,M) = (Tr [S'! (U - ^) 5 (U + ^)]"). Ex¬ 
panding the trace, we find a multiple sum over trajectories. 


Cn{e,M) 


1 ” 

n y y ax 

A/fl^n 1 1 Z—/ Z—/ I ^ 

H fc=l 2fc,Ofc 7/c,crfc 


So 


2^(U+ra) 




(9) 


such that 7 fc goes from to Ok, while goes from to Ofc+i- The channels labels are all 
being summed from 1 to M. 

In (|9]) we have used 

fh fh 

( 10 ) 

where Ty is the total duration of 7 . The quantity = Ofc A^^ is a collective stability, while 
y = T,kAk A ^ Y,kAk collective action and duration of the 7 trajectories, and 

analogously for a. 

The result of the sum in ([9]) is, for a chaotic system, a strongly fluctuating function of the 
energy. The average over E, under the stationary phase approximation, requires 7 and a to 
have almost the same collective action. In the past years^, it has been established that these 
action correlations arise when each a follows closely a certain 7 for a period of time, and 
some of them exchange partners at so-called encounters. A g-encounter is a region where q 
pieces of trajectories run nearly parallel and q partners are exchanged. This theory has been 
presented in detail before^^i^ We consider only systems not invariant under time-reversal, 
so a trajectories never run in the opposite sense with respect to 7 trajectories. 

For example, we show in Figure la a situation contributing to the second correlation 
function, C 2 {e, M). Trajectory 71 starts in channel ii and ends in channel Oi, while 72 starts 
in channel *2 and ends in channel 02 . On the other hand, cJi and cr 2 are initially almost 
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FIG. 1 . a) Correlated trajectories contributing to C2{e,M). Solid lines are 71 (going from ii to 
oi) and 72 (going from 12 to 02), dashed lines are cri (going from ii to 02) and cj 2 (going from 12 
to oi). In this situation we have one 2 -encounter and one 3 -encounter (the encounters are greatly 
magnified), b) Correlated trajectories contributing to C3(e,M), in a case with coinciding channels. 
In both figures the chaotic nature of the trajectories is not shown. 


identical to 71 and 72, respectively, but they exchange partners in a 2 -encounter. Later, 72 
has a 3 -encounter with itself, inside of which the pieces of cti are connected differently. We 
also show in Figure lb a situation contributing to 6*3(6, M) which has no encounters, but has 
coinciding channels. There are two major simplihcations done here for visual clarity: 1 ) The 
encounters are greatly magnihed, to show their internal structure; 2 ) The actual trajectories 
are extremely convoluted and chaotic. Many other examples of correlated trajectories can 
be found in previous work^^i2,“Jiiii£i^ 

Correlated sets of trajectories contributing to the semiclassical calculation of correlation 
functions can be depicted in the form of ribbon graphs, as suggested in^i^. The g-encounters 
become vertices of valence 2 q. Channels also become vertices, but their valence depends on 
whether there are coinciding channels or not. The pieces of trajectories connecting vertices 
become fat edges, or ribbons. Each ribbon is bordered by one 7 and one a, and these 
trajectories traverse the encounter vertices in a well dehned rotation sense: a trajectory 
arriving from one ribbon departs via the adjacent ribbon (graphs endowed with a cyclic order 
around vertices are also called maps). We show in Figure 2 the ribbon graphs corresponding 
to the trajectories shown in Figure 1 . 

Following previous work on transport and on closed systems, Kuipers and Sieber obtained 
some diagrammatic rulesi^, that determine how much a given graph contributes to the cor¬ 
relation function. The contribution of a graph factorizes into the contributions of individual 
vertices and edges: an encounter vertex of valence 2 q gives rise to -M(l - iqe)] channels of 
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FIG. 2. The ribbon graphs corresponding to Figure 1. Each ribbon is bordered by one 7 and one 
a. Ribbons only meet at vertices, and g-encounters become vertices of valence ‘Iq. 


any valence give rise to M; each ribbon gives rise to [M(l - These rnles were then 

used in several works dealing with time delay statistics^ii^iSZiS^. 

Notice that there are no periodic orbits in a ribbon graph that arises from the semiclassical 
expansion of time delay. This means that we may start from and follow cxi up to oi, then 
follow 7 i in reverse back to ii, then (T 2 to 02 , then 72 in reverse back to Z 2 , and so on, and 
traverse every border of every ribbon exactly once. This means that the graph has a single 
face. 

The contribution of a graph will be proportional to where V is the total number 

of vertices (including channels) and E is the total number of edges. The Euler characteristic 
of a ribbon graph is V - E + F, where F is the number of faces (E = 1 in our case). The 
Euler characteristic is also equal to 2 - 2g, where g is called the genus. Therefore, the 
1/M expansion coming from semiclassical diagrammatics is actually what is called a genus 
expansion: the contribution of a graph is proportional to l/M^^. Graphs with g = 0 are 
called planar (they can be drawn on the plane so that the ribbons never cross each other), 
and they give the leading order contribution. 

The graph in Figure 2a, for example, contributes 

(l-2fe)(l-3if) 

to ( 72 ( 6 , M). Notice that it is not a planar graph, since there is a crossing between two of 


7 




















the ribbons. This particular graph actually has g = 1 (this means it may be drawn on a 
torus without any crossings). The graph in Figure 2b, on the other hand, is planar and 
contributes to C 3 {e,M). 


B. Gaussian integrals and Wick diagrammatics 


We shall introduce a certain Gaussian matrix integral and formulate it diagrammatically, 
using Wick’s rule. This procedure has been discussed in detail for hermitian matricesJ^i^ 
The only difference compared to the present work is that we integrate over non-hermitian 
matrices. Our diagrams are then interpreted as providing the semiclassical formulation of 
the time delay problem. The same approach has been used to treat transport statisticsj^Si 
Let Z denote a general complex matrix of dimension N, and dehne 

where the normalization constant is 

Z = J (13) 

We see (1121) as an average value, but we use the symbol ((■)) to differentiate it from the true 
physical energy-average we considered in previous sections. For example, since the elements 
are actually independent, it is clear that 


Integrals over a product of matrix elements can be computed using the so-called Wick’s 
rule, which states that we must sum, over all possible pairings between Z’s and Zl’s, the 
product of the average values of the pairs. Namely, 


n 

\A:=1 


^mkjk^lkrk , 


n 

a^Sn fc=l 


(15) 


If we the quantity we wish to average involves traces of ZZl, all we need to do is expand 
these traces in terms of matrix elements and apply Wick’s rule. Most importantly, we can 
then employ a diagrammatic technique. 

For example, suppose we wish to compute 


((Tr[(ZZt)2]Tr[(zzt)3]^_____ 


(16) 
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We start by writing it as 

E E (17) 

mi,...,m5 W U=1 s=3 J II 

where all sums run from 1 to (in the hrst product we mean m 3 = mi, while in the second 
product we mean mg = m 3 ). The diagrammatics consists in picturing the matrix elements as 
pairs of arrows. Arrows that represent elements from Z have a marked end at the head, while 
arrows that represent elements from Zt have a marked end at the tail. Arrows representing 
matrix elements coming from traces are arranged in clockwise order around vertices, so that 
all marked ends are on the outside. Finally, the elements that do not come from traces 
are arranged surrounding the other ones, also in clockwise order. Since this is most easily 
explained by means of an image, we show it in Figure 3(a). 

Once we have arranged the arrows, Wick’s rule consists in making all possible connections 
between them, using the marked ends. Clearly, this produces a ribbon graph. According 
to Eq. ffTTj) . when computing the value of a graph, each ribbon gives rise to a factor 
For the example in Figure 3(a), there are 7! possible connections. We show two of them in 
Figures 3(b,c). The coupling in Figure 3(b) leads to the identihcations 

4= mi, i2= 7712 =1713 =1714 = 1715 Oi = j2 = js = ji = jb, O 2 = jl, (18) 

and gives a contribution of to the average flTbl) . Notice how this coupling is similar to 

Figure 2. On the other hand, the coupling in Figure 3(c) leads to the identihcations 

il=7ni, ^2 = m 2 =7714 = 7715 01=^2=23=75, O 2 = jl- (19) 

In this case the indices m 3 and remain free to be summed over. Therefore, this coupling 
gives a contribution of to the average flTHl) . 

Free indices arise from closed loops in the ribbon graph. Each such loop increases by one 
the number of faces of the graph (every graph has at least one face). Therefore, the power 
of N in the contribution of a given coupling is always one less than the number of faces in 
the graph. 

It should be clear that this theory is very close to the semiclassical approach to time 
delay, provided we choose hi = M(1 -ie). However, the ribbon graphs in the semiclassical 
theory always have a single face. As we have just mentioned, this corresponds to keeping 
only those Wick couplings whose contribution does not depend on N. Since all contributions 
are proportional to a positive power of N, we can simply let A^ ^ 0. 
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FIG. 3. Diagrammatics of Wick’s rule, for the average in (|16ll . In a) we see how the matrix 
elements are turned into arrows with marked ends and, in the case of traces, arranged clockwise 
around vertices. In the vertex of valence 6 we have written each label only once, for clarity. In b) 
and c) we see two particular Wick couplings, out of the possible 7!. The labels of the arrows in b) 
and c) are the same as in a). Notice the similarity between b) and Figure 2. 


C. Matrix integrals for correlation functions 


Let ^ = (12---n) be the cyclic permutation of the hrst n positive integers, and let i = 
(zi,...,z„) and o= (oi,...,o„). Introduce the integral 

1 






( 20 ) 


This can be seen as a Gaussian average as the ones considered previously, if we understand 
the hrst term in the exponent, measure. Accordingly, we 

set 

Z = J (21) 

The rest of the exponential can be Taylor expanded as 

_ £ {-My (1 - Tr[(ZZt)g] 

t=0 \g>2 Q 



10 


































For now, we consider this as a formal power series and integrate term by term, employing 
Wick’s rule and its diagrammatical representation previously discussed. By construction, 
encounter vertices of valence 2q will be accompanied by the factor -M(l - iqe), giving the 
correct semiclassical diagrammatic rules. 

The integral fl20|) is therefore designed to automatically produce all the required ribbon 
graphs for the semiclassical evaluation of the correlation function Cn{^, M). The exponential 
produces all possible encounters, while the matrix elements in the last product play the role 
of the channels. In line with Eq.([n]), we must sum over all channels from 1 to M, i.e. we 
must consider the quantity 

M M 

gn{e,M,N) = J^Gnie,M,N,ld)= ^ ^ G^ie, M, N,ld). (23) 

2jO ^ 1 5’'T^n = l = l 

The matrix integral produces more graphs than needed, but we have provided for this 
overcounting. For example, the Taylor series of the exponential naturally has a t! in the 
denominator, which is responsible for eliminating the symmetry associated with shuffling 
the vertices, when there are t of them. Also, graphs are produced that differ from each other 
only by the rotation of a vertex. This is why we have divided Tr[(ZZl)'J] by q: it remedies 
the overcounting that would be caused by the possible q rotations of the vertex. 

As we have discussed, in order to select only those ribbon graphs with a single face it is 
necessary to take the limit A^ 0 at the end of the calculation. Therefore, the correlation 
function will be given by 

a(e,M) = lim61„(e, M,iV). (24) 

It is not very difficult to implement Eq. fl20l) in a computer and obtain the hrst few orders 
in 1/M for the hrst few correlation functions (the integral is not to be performed numerically, 
of course, but using Wick’s rule together with the covariance (ITT)) h This leads to the results 
in ([S])-(I7j). Notice that letting ^ 0 in this context presents no difficulty. 


IV. EXACT SOLUTION 

This Section is dedicated to the exact solution of the matrix integral fl20|) . and the cal¬ 
culation of its limit as N ^ 0. 
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A. Angular integration 


Introduce the singular value decomposition Z = UDV, where D is real, positive and 
diagonal while U and V are unitary. Let X = be a matrix with the same eigenvalues as 
ZZ\ and denote these eigenvalues by Xj, 1 < i < N. It is known^ that the measure dZ is 
expressed in these new variables as 

dZ = CN\A{X)\‘^dxdUdV, (25) 

where cn depends only on the dimension, dU is the normalized Haar measure on the unitary 
group h{{N), and the Vandermonde squared is the Jacobian of the transformation. This is 
a generalization of the transformation from cartesian to polar coordinates in the complex 
plane. We shall hrst perform the angular integration over U and V. 

A minor point to be mentioned is that dV is not the same as the normalized Haar 
measure. This is related to the fact that in the singular value decomposition there is a 
certain ambiguity, as we may freely conjugate D by a diagonal unitary matrix. The matrix 
V is thus uniquely determined only as an element of the coset U{N)/[U{1)]^. However, the 
functions we shall integrate, polynomials in matrix elements as those in Section 3.3, are all 
invariant under multiplication by a diagonal unitary matrix, and in this context dV behaves 
just like the Haar measure, up to normalization. 

The only part of the integral in fl20|) that depends on the angular variables U and V is 
the last product. Thus, the angular integral we need is 

/ n 

dUdVW y (26) 

k=ljk,mk 

Given j = ■■■jn), m = {mi,m 2 , ...,mn) and r e Sn, let 

n 

k=l 

The Weingarten function of the unitary group is dehned by 

/ n n 

<1UUU.aUU= E (28) 

k=l a,T^Sn k=l 

and its character expansion is22i~— 

Wg^(^) = ^ E J^xx{g). (29) 

£{\)<N 
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Using the above machinery, we have 


A= X! ^)Wg^(rcr ^)p^-ig{X)S^[i,i]6p[o,^{o)], (30) 

arp6£Sn 

where px are the power sum symmetric functions, and we have used that 

n n 

n E Dj^D^^5r[3,rn\5B[j,m] = Y\Y,^3A-^e[3.j]=Pr-^e{X). (31) 

k=ljk,mk fc=l jk 

The quantity we are after, Eq. fl2^ . requires summation over the indices i and o. It is 
easy to see that 

M 

£ 5.[z,z]=M^('^)=p.(l^), (32) 

where i{a) denotes the number of cycles of the permutation a, and 

M 

£ 5,[o,e(o)] = M^(^«)=p,^(l^). (33) 

01,--,On = l 


Notice that the channel labels in the original matrix integral fl2U]) are all constrained to be 
between 1 and N. Nevertheless, we are summing them from 1 to M. We are thus assuming 
N > M. However, this will not deter us from letting ^ 0 later. 

Once we expand 

p^MX)=Y.Xx{r-^0)sx{X), (34) 

Aun 

where sx are the Schur symmetric functions, we get 

E^=E E Wg^(pr^)Wg^(Tcr^^)XA(W^0)sA(A:)p<,(l^)pp^(l^). (35) 

1 5 X\-n arpO^Sn 

Repeated use of the character orthogonality relation 

___ 77 1 

E Xik(^)X\{ra) = —XA(a)(5^,A (36) 

reSn «A 

leads to ^ 

B. Eigenvalue integration 

So far, the quantity we are after is given by 

5„(€, M, Af) = y G„ = y XaK) j K(e, M, A^), (38) 


13 




where Tl{e, M, N) is the radial integral over the eigenvalues of ZZ^. It is equal to 

n{e,M,N) = ^ rd5det(l-X)^e"*^T’'[T^]|A(x)psA(X), 

MZ Jo 

where we have used that 


-ME„. ^ 


From the well known Schur function expansion, 

(Mie) 


m=0 /ii-m ^ ^ 


we get 


where 


°° ( 

-A E d,I,JM,N), 


m=0 


m\ 


fi\-m 


= dxdetil - X)^\A{x)\h, s,{X). 

Using the determinantal form of the Schur functions and the identity 

A| 


(i^) 


A(A) 


one can show that 


^A,/. - 


Cat A! 


det 


det(l-X)^-i’ 

(M - +j-l)\{Xi-i + fij -j + 2N)\' 


{M + 2N + \i-i)\ 

Two factorials can be taken out of the determinant, and we can write 

c nN\ ^ {M-iXj+j - 1 )! 

Z l{{M + 2N + X,-j)\ 


1x,ti - 


— det ((Ai -i + fXj-j + 2N)\). 


Introducing (M + j - 1)! in the product, we get 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 


C. The A ^ 0 limit 


We must now take the A ^ 0 limit. This is a delicate procedure. We can only do 
it for quantities that are analytic functions of A. For example, using the singular value 
decomposition, the normalization constant fl?I|) becomes 

r' oo p N 

^ = cnJ^ -ypjnj!(A-j)!. (48) 
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It is perfectly fine to take the limit in the denominator. In the rest of the expression, we 
must leave N intact for now. In this sense, we write 


N 




(49) 


The quantity contains the factor 


^ {M + N-j)\ 

|.l(M + 2Af + A,-i)! 


(50) 


First, we let 0 inside the product, to get 




(51) 


This still depends on N via the limit of the product. However, Xj = 0 for j > £(X). Hence, if 
we assume N > £(X), we can write this as 


'<A) (M-j)! 


(62) 


P(M + A,-j)! [MY' 

which is independent of N. Now, in all rigor we are not allowed to take iV ^ 0 after assuming 
N > i{\). We do it anyway, and write 




N\ 


■det ((Ai -i + j + 2N)\). 


(53) 


Z [M]^[MY 

Further, we factor out the smallest factor from each row of the determinant, producing 


nj=i(^ + - j + Pat)!. If we assume that N > £{fi), then = 0. Hence, using 

m « (jVtA,-j)! / (A.-^^-27V)! \ 

" [M],.[Af]yU {N-j)\j\ ““1 (A,-i + ]V)! j 

We again consider N > £{X) hrst and ^ 0 later, to arrive at 


(54) 


I 


A,/i 


[Nf 


[MY[MY 


det 


{\i - i + fij - j + 2N)\ 

, (Ai-i + iV)! , 


(55) 


1. The determinant 


We need to consider the determinant 




( 56 ) 


16 














where 


(57) 


ai = \i - i + N, bj = fj,j - j + N. 

Each column consists of raising factorials, i.e. we have 

(a* + l)(ai + 2)---(aj + bj) = -—^ -. (58) 

di 

We therefore expand each column using identity 


[U" = E 

k=0 


X^, 


in terms of unsigned Stirling numbers of the first kind. We get 




a. 


bi + l bj+l 

= E 

k-i=l 


bj + 1 

II 

1 

bj + 1 

kj 

kj=0 

kj + 1 


k-i 

a/. 


The determinant is then given by 


N bj 

»=nE 

j=l kj=0 


bj + 1 
kj +1 


det(af^). 


Introducing kj = Uj - j + N we have 


(59) 


(60) 


(61) 


N H 

p-n E 


j=l uij=j-N 


^j-j + N + l 
Uj — j + N + \ 




(62) 


Notice that u is not a partition, since its elements are not necessarily ordered, and they 
can be negative. Still, the last determinant, if it does not vanish, can be turned into a Schur 
function by simply re-ordering the columns. Let u be the partition that is created in this 
way, and \u\ the number it partitions. For instance, if a; = (1,1,-1,1) we have 

det ^af af“^ af"^ af~^ af"®-••) = -det ^af af~^ ’ (63) 


so the corresponding partition is a; = (1,1) and |c(;| = 2. As we can see, the reordering of the 
columns may lead to a change in sign. Let rji^u) denote this sign, so that 

det = r](u)A(a)sz{a) = r 7 ((:u)^[A^]^ss(a) Yl J- (64) 

^ ' nl j^\ 

We must consider the 0 limit of 

sM = E \Cp\xMPp(a). (65) 

1 ^ 1 - 


16 



Fx,p 


h 


(1) 

1 

1 

1 

2 

1 

2 

6 

2 

0 

2 

6 


(11) 

2 

4 

2 

12 

4 

4 

48 

12 

0 

8 

12 


(2) 

2 

2 

4 

4 

4 

12 

12 

8 

0 

12 

48 


(111) 

3 

9 

3 

36 

9 

6 

180 

36 

0 

18 

18 


(21) 

3 

6 

6 

18 

12 

18 

72 

36 

0 

36 

72 

A 

( 3 ) 

3 

3 

9 

6 

9 

36 

18 

18 

0 

36 

180 


(1111) 

4 

16 

4 

80 

16 

8 

480 

80 

0 

32 

24 


(211) 

4 

12 

8 

48 

24 

24 

240 

96 

0 

72 

96 


(22) 

4 

10 

10 

32 

28 

32 

132 

92 

12 

92 

132 


(31) 

4 

8 

12 

24 

24 

48 

96 

72 

0 

96 

240 


( 4 ) 

4 

4 

16 

8 

16 

80 

24 

32 

0 

80 

480 


TABLE I. The function for the first few values of A and fi. 


The limit of Pp{a) can be obtained simply removing from this quantity everything that scales 
with N: 


i(p) 

hmpp{{Xi-i +N}) = n 


e(\) 

i=l 



-/p(A). 


( 66 ) 


We can hnally write 


V 



where the function Fx^p is given by 


dx 


n\ 




(67) 


F: 


A,/i. 


n E 


j=l U)j=j-l{p.) 


hi - i +1 
Uj -j + l 


1 ^1- p^\uj\ 


( 68 ) 


We tabulate some values of this function in Table 1. Several properties stand out from 
inspection of these values. First, all values are non-negative integers. Second, and this is 
easy to prove, that 

^A,(i) =/(i)(A) = n if A K n. (69) 


Third, and this we can only conjecture, that 




( 70 ) 
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where X' and jj.' are the conjugate partitions of A and fi, respectively. Finally, we notice 
that, when p is not a hook partition, it seems that is different from zero only if A is 
also not a hook partition. We have checked this extensively, but do not have a proof. 


D. Final Result 


It is time to put the pieces back together. We have to plug the limiting value of V into 
the expression for Eq. (j55|l . put this into the expression for the radial integral, Eq. (02]), 
and hnally arrive at the quantity we want, which is On, Eq. fIdSp . After some cancelations, 
we get that the limit as N ^ 0 of On, which is nothing but the semi classical expression for 
the correlation function Cn{^,M), is given by 

[MY 


1 °° (MifY 

\na g,Xe-,M,N) = C„(e, A/) = — ^ 1 

N-*o Mn\ ml 


fi\-m X\-n 




(71) 


This expression is perhaps not as simple we one might hope for, specially the Fa_^ part. 
This complication is probably due to the fact that we are using a Taylor series in e. We know 
that, at each order in 1/M, the correlation functions are rational functions of e, with the 
denominator being a power of (1 - ie). Maybe if this fact could be explicitly incorporated 
into the calculation somehow, the resulting expression would be more manageable. 

We can see from (16^ that Cn{e, M) = 1 + me + O(e^). 

Using ([2|), our semiclassical approach leads to the following expression for the average 


value of time delay moments Aim = igTr/Q™): 


{Mm) - E 


di. 


ml 


fi\-m 


[M] 


E 


(- 1 )' 


M n=l 


n\ 


(’") E (72) 

V ^ / Xi-n 


On the other hand, we have shown^ using the random matrix theory approach, that 

di, 


{Mm) 


RMT - 


ml 


E 

/ii-m 


[Ml 




(73) 


Therefore, if the identity 

"7 Xm-n /™\ 

E ^ „ ) E dxXx{n)[MYFx,f, = Xt^{m)[MY (74) 

holds, the semiclassical formula for {Mm) becomes exactly equal to the corresponding RMT 
prediction. We have checked that fl74|) indeed holds for all p f- m up to m = 8. This 
guarantees agreement between the semiclassical and RMT calculations up to the first 8 
moments. We expect that this agreement should in fact hold for all moments, and in fact 
for all polynomial functions of the time delay matrix. 
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